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Abstract. Non radial oscillations of neutron stars are associated with the emission of 
gravitational waves. The characteristic frequencies of these oscillations can be com- 
puted using the theory of stellar perturbations, and they are shown to carry detailed 
information on the internal structure of the emitting source. Moreover, they appear to 
be encoded in various radiative processes, as for instance in the tail of the giant flares 
of Soft Gamma Repeaters. Thus, their determination is central to the theory of stellar 
perturbation. A viable approach to the problem consists in formulating this theory as 
a problem of resonant scattering of gravitational waves incident on the potential bar- 
rier generated by the spacetime curvature. This approach discloses some unexpected 
correspondences between the theory of stellar perturbations and the theory of quantum 
mechanics, and allows us to predict new relativistic eff"ects. 

Keywords : gravitational waves - black hole physics - stars: oscillations - stars: neu- 
tron - stars: rotation 



1. Introduction 

The theory of stellar perturbations is a very powerful tool to investigate the features of gravi- 
tational signals emitted when a star is set in non-radial oscillations by any external or internal 
cause. The characteristic frequencies at which waves are emitted are of great interest in these 
days, since gravitational wave detectors Virgo and LIGO are approaching the sensitivity needed 
to detect gravitational waves emitted by pulsating stars. These frequencies carry information on 
the internal structure of a star, and appear to be encoded in various radiative processes; thus, their 
study is a central problem in perturbation theory and in astrophysics. In this paper I will illus- 
trate the theory of perturbations of a non-rotating star, describing in particular the formulation 
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that Chandrasekhar and I developed, its motivation and outcomes. Furthermore, I will briefly 
describe how the theory has been applied to study the oscillation frequencies of neutron stars. 

In order to frame the problem in an appropriate historical perspective, it is instructive to re- 
mind ourselves how the study of stellar perturbations was treated in the framework of Newtonian 
gravity. In that case, the adiabatic perturbations of a spherical star are described by a fourth-order, 
Unear, diff'erential system which couples the perturbation of the Newtonian potential to those of 
the stellar fluid. All perturbed quantities are Fourier-expanded and, after a suitable expansion 
in spherical harmonics, which allows for the separation of variables, the relevant equations are 
manipulated in such a way that the quantity which is singled out to describe the perturbed star 
is the Lagrangian displacement ^ experienced by a generic fluid element; indeed, the changes 
in density, pressure and gravitational potential induced by the perturbation, can all be expressed 
uniquely in terms of ^. The equations for ^ have to be solved by imposing appropriate boundary 
conditions at the centre of the star, where all physical quantities must be regular, and on its bound- 
ary, where the perturbation of the pressure must vanish. These conditions are satisfied only for a 
discrete set of real values of the frequency, {aj„), which are the frequencies of the star's normal 
modes. Thus, the linearized version of the Poisson and of the hydro equations are reduced to a 
characteristic value problem for the frequency o). 

An adequate base for a rigorous treatment of stellar pulsations of a spherical star in general 
relativity was provided by K.S. Thorne and collaborators in a series of papers published in the 
late sixties-early seventies of the last century (Thorne & Campolattaro 1967, 1968; Campolattaro 
& Thorne 1970; Thome 1969a,b; Ipser & Thome 1973). The theory was developed in analogy 
with the Newtonian approach, and was later completed by Lindblom & Detweiler (1983), who 
brought the analytic framework to a form suitable for the numerical integration of the equations, 
thus allowing for the determination of the real and imaginary part of the characteristic frequencies 
of the ( - 2, quasi-normal modes. Indeed, a main difference between the Newtonian and the 
relativistic theory is that in general relativity the oscillations are damped by the emission of 
gravitational waves, and consequently the mode eigenfrequencies are complex. Higher order 
{€ > 2) mode frequencies were subsequently computed by Cutler & Lindblom (1987). 

In 1990, Professor Chandrasekhar and 1 started to work on stellar perturbations, and we 
decided to derive ab initio the equations of stellar perturbations following a difi'erent approach, 
having as a guide the theory of black hole perturbations rather than the Newtonian theory of 
stellar perturbations. 



1.1 Black hole perturbations: wave equations and conservation laws 

In 1957 T. Regge and J. A. Wheeler set the basis of the theory black hole perturbations showing 
that, by expanding the metric perturbations of a Schwarzschild black hole in tensorial spherical 
harmonics, Einstein's equations can be separated (Regge & Wheeler 1957). Spherical harmonics 
belong to two dififerent classes, depending on the way they transform under the parity transfor- 
mation — » TT - and ^ — » TT -H ^; those which transform like (-1)'^'^+^' are named odd, or axial. 
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those that transform like (-1)'^ are named even, or polar. The perturbed equations decouple in 
two distinct sets belonging to the two parities. Regge & Wheeler further showed that by Fourier- 
transforming the time dependent variables, the equations describing the radial part of the axial 
perturbations can easily be reduced to a single Schroedinger-like equation, and 13 years later 
Frank Zerilli showed that this can also be done for the much more complicated set of polar equa- 
tions (Zerilli 1970a,b); thus, the axial and polar perturbations of a Schwarzschild black hole are 
described by the wave equation 
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wherer, = rH-2Mlog^^^ - lj,n = i(^-Hl)(/'-2), and M is the blackhole mass. The superscript 

- and + indicate, respectively, the Regge- Wheeler equation for the axial perturbations, and the 
Zerilli equation for the polar perturbations, and the corre spending potentials. The Regge- WTieeler 
potential for ^ = 2 is shown in figure [1] 




Figure 1. The potential barrier generated by the axial perturbations of a Schwarzschild black hole for 1 = 2. 
The potential for the polar perturbations has a similar form. 



An alternative approach to studying black hole perturbations considers the perturbations of 
the Weyl and Maxwell scalars within the Newman-Penrose formalism. Using this approach in 
1972 S. Teukolsky was able to decouple and separate the equations governing the perturbations 
of a Kerr black hole (Teukolsky 1972, 1973), and to reduce them to a single master equation for 
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the radial part of the perturbation Ri,„: 



ARlm,rr + 2(s + l)(r - M)Ri,n,r + V{oj, r)Ri,„ = , 



2Mr + c?-. 



(4) 



Variable separation was achieved in terms of oblate spheroidal harmonics, and the potential 
r) is given by 



where a is the black hole angular momentum, A;,„ is a separation constant, and s, the spin-weight 
parameter, takes the values s - 0, +1,±2, respectively for scalar, electromagnetic and gravita- 
tional perturbations. It may be noted that for the Schwarzschild perturbations, due to the back- 
ground spherical symmetry, non-axisymmetric modes with a e""* dependence can be deduced 
from axisymmetric, m - modes by suitable rotations of the polar axes. As a consequence, 
the potentials (|2|l and (O do not depend on the harmonic index m. Conversely, since Kerr's 
background is axisymmetric, this degeneracy is removed and the potential (|6]l depends on m. 
Moreover, while the Schwarzschild potentials (|2) and (l3) are real and independent of frequency, 
the potential barrier of a Kerr black hole is complex, and depends on frequency. 

The wave equations governing black hole perturbations show that the curvature generated 
by a black hole appears in the perturbed equations as a one-dimensional potential barrier; con- 
sequently, the response of a black hole to a generic perturbation can be studied by investigating 
the manner in which a gravitational wave incident on that barrier is transmitted, absorbed and 
reflected. Thus, the theory of black hole perturbations can be formulated as a scattering the- 
ory, and the methods traditionally applied in quantum mechanics to investigate the behaviour 
of physical systems described by a Schroedinger equation can be adapted and used to study the 
behaviour of perturbed black holes. For instance, it is known that in quantum mechanics, given 
a one-dimensional potential barrier associated with a Schroedinger equation, the singularities in 
the scattering cross-section correspond to complex eigenvalues of the energy and to the so-called 
quasi-stationary states. Since the perturbations of a Schwarzschild black hole are described by the 
Schroedinger-like equation ([U with the one-dimensional potential barriers (|2} and (l3), in which 
the energy is replaced by the frequency, the singularities in the scattering cross-section will pro- 
vide the complex values of the black hole eigenfrequencies, and the corresponding eigenstates 
will be the black hole Quasi Nomal Modes (QNM). These modes satisfy the boundary condi- 
tions of a pure outgoing gravitational wave emerging at radial infinity (r» — > -Hoo), and a pure 
ingoing wave impinging at the black hole horizon (r* — > -oo). That these solutions should exist 
had been suggested by C.V. Vishveshwara in 1970 (Vishveshwara 1970), and the next year W.H. 
Press confirmed this idea by numerically integrating the axial wave equation ([T), and by showing 
that an arbitrary initial perturbation ends in a ringing tail, which indicates that black holes pos- 
sess some proper modes of vibration (Press 1971). However, it was only in 1975 that S. Chan- 
drasekhar and S. Detweiler computed the complex eigenfrequencies of the quasi-normal modes 
of a Schwarzschild black hole, by integrating the Riccati equation associated with the axial equa- 
tion ([T]| (Chandrasekhar & Detweiler 1975). In addition, they also showed that the transmission 
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and the reflection coefficients associated respectively with the polar and with the axial potential 
barriers are equal. As a consequence, the polar and the axial perturbations are isospectral, i.e. the 
polar and axial QNM eigenfrequencies are equal. This equality can be explained in terms of a 
transformation theory which clarifies the relations that exist between potential barriers admitting 
the same reflection and absorption coefficients. This is an example of how the scattering approach 
has been effective not only to determine the QNM frequencies, but also to investigate the inner 
relations existing among the axial and polar potential barriers and to gain a deeper insight in the 
mathematical theory of black holes, which was illustrated by S. Chandrasekhar in his book on 
the subject (Chandrasekhar 1984). Following this approach, a variety of methods developed in 
the context of quantum mechanics have been used to determine the QNM spectra of rotating and 
non-rotating black holes, such as the WKB and higher order WKB method, phase-integral meth- 
ods and the theory of Regge poles, just to mention some of them (Schutz & Will 1985; Ferrari 
& Mashhoon 1984a,b; Andersson, Araujo & Schutz 1993a,b,c; Andersson 1994; Andersson & 
Thylwe 1994). 



1.2 A conservation law for black hole perturbations and its generalization to perturbed 
stars 

In quantum mechanics the equation 

|/?|2 + |r|2 = 1, (7) 

where R and T are the reflection and transmission coefficients associated with a potential barrier, 
expresses the symmetry and unitarity of the scattering matrix; it says that, if a wave of unitary am- 
plitude is incident on one side of the potential barrier, it gives rise to a reflected and a transmitted 
wave such that the sum of the square of their amplitudes is still one. Therefore, equation ^ is an 
energy conservation law for the scattering problem described by the Schroedinger equation with 
a potential barrier This conservation law is a consequence of the constancy of the Wronskian 
of pairs of independent solutions of the Schroedinger equation. Similarly, the constancy of the 
Wronskian of two independent solutions of the black hole wave equations allows us to write the 
same relation between the reflection and transmission coefficients associated with the potential 
barrier, and therefore it shows that such an energy conservation law also governs the scattering of 
gravitational waves by a perturbed black hole. It should be stressed that such energy conservation 
law does not exist in the framework of the exact non-linear theory; however, it can be derived in 
perturbation theory both for Schwarzschild, Kerr and Reissner-Nordstrom black holes. 

This possibility led Chandrasekhar to the following consideration. Since in general relativity, 
any distribution of matter (or more generally energy of any sort) induces a curvature of the space- 
time - a potential well - instead of picturing the non-radial oscillations of a star as caused by 
some unspecified external perturbation, we can picture them as excited by incident gravitational 
radiation. Viewed in this manner, the reflection and absorption of incident gravitational waves 
by black holes and the non-radial oscillations of stars, become different aspects of the same ba- 
sic theory. However, this idea needed to be substantiated by facts, and our starting point was 
to show that also for perturbed stars it is possible to write an energy conservation law in terms 
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of Wronskians of independent solutions of the perturbation equations of a spherical star. This 
is easy if we consider the axial perturbations of a non-rotating star, because in that case, as we 
shall show in Section |2] the perturbed equations can be reduced to a wave equation with a one 
dimensional potential barrier as for black holes. However, to derive the conservation law for the 
polar perturbations was not easy, because the corresponding equations are a fourth order linear 
differential system, in which the perturbed metric functions couple to the fluid perturbations, and 
it was not clear how to define the conserved current. Anyway, working hard on the equations, 
we were able to derive a vector ^ in terms of metric and fluid perturbations, which satisfies the 
following equation (Chandrasekhar & Ferrari 1990a): 

—E" = 0, a^(x^ = r, = §). (8) 

(It is worth reminding ourselves that, due to the spherical symmetry, it is not restrictive to con- 
sider axisymmetric perturbations m - 0). The vanishing of the ordinary divergence implies that, 
by Gauss's theorem, the flux of across a closed surface surrounding the star is a constant. 
When the fluid variables are switched off, this conservation law reduces to that derived for a 
Schwarzschild black hole, and therefore we thought we were on the right track. However, there 
was still a question to answer: are we entitled to say that the vector E actually represents the flux 
of gravitational energy which develops through the stars and propagates outside? If so, equation 
^ should reduce to the second variation of the time component of the well known equation 

A[V=i(r'^^ + n]=0. (9) 

where f is the stress-energy pseudotensor of the gravitational field. The problem is that f is 
not uniquely defined; indeed equation (|9]i shows that it is defined up to a divergenceless term. 
A possible definition is that given by Landau & Lifschitz (1975), which has the advantage of 
being symmetric. However, the second variation of the time component of equation (|9]l assum- 
ing f - does not give the divergenceless equation satisfied by our vector E, neither for 
the Einstein-Maxwell case, nor in the case of a star Then, Raphael Sorkin suggested that the 
pseudo-tensor whose second variation should reproduce our conserved current is the Einstein 
pseudo-tensor, because its second variation retains its divergence-free property, provided only 
the equations governing the static spacetime and its linear perturbations are satisfiec(3. This prop- 
erty is a consequence of the Einstein pseudo-tensor being a Noether operator for the gravitational 
field; the Landau-Lifshitz pseudotensor failed to reproduce the conserved current because it does 
not satisfy the foregoing requirements. In addition, Sorkin pointed out that the contribution of 
the source should be introduced not by adding the second variation of the source stress-energy 
tensor T'^'', as one might naively have thought, but through a suitably defined Noether operator, 
whose form he derived for an electromagnetic field (Sorkin 1991). Though this operator does not 
coincide with T'^'', it gives the same conserved quantities. Thus, the flux integral which we had 
obtained, I would say, by brute force, working directly on the perturbed hydrodynamical equa- 
tions, could be obtained from a suitable expansion of the Einstein pseudo-tensor showing that. 



It should be mentioned that the first variation of the Einstein pseudo-tensor vanishes identically. 
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as for black holes, energy conservation also governs phenomena involving gravitational waves 
emitted by perturbed stars (Chandrasekhar & Ferrari 1991a). We therefore decided to derive ab 
initio the equations of perturbations of a spherical star in the same gauge used when studying the 
perturbations of a Schwarzschild black hole, and to study the problem as a scattering problem. 
In the next sections I shall briefly illustrate the main results we obtained by using this approach 
(Chandrasekhar & Ferrari 1990b, 1991b,c, 1992; Chandrasekhar, Ferrari & Winston 1991). 



2. Perturbations of a non-rotating star 

As for a Schwarzschild black hole, when the equations describing the perturbations of a spheri- 
cal star are perturbed and expanded in spherical tensor harmonics they decouple in two distinct 
sets, one for the polar and one for the axial perturbations. The axial equations do not involve 
fluid motion except for a stationary rotation, while the polar equations couple fluid and metric 
perturbations. The axial equations are therefore much simpler and we showed that, after separat- 
ing the variables and Fourier-expanding the perturbed functions, they can be combined as in the 
Schwarzschild case, and reduced to a single Schroedinger-like equation with a one-dimensional 
potential barrier (Chandrasekhar & Ferrari 1990b, 1991c): 

_A + [^2 _ VJirWJ = 0, (10) 

where ^ 

r. = r e-^^f'-dr, (11) 
Jo 

and 



VJ(r) = —[{({ + l)r + r\€ - p) - 6m(r)]. (12) 
r 

The functions v(r) and yU2(r), which appear in the definition of the radial variable r,, are two 
metric functions which are found by solving the equations of stellar structure for an assigned 
equation of state (EOS). e(r) and p{r) are the energy density and the pressure in the unperturbed 
star; outside the star they vanish and equation (fT2l i reduces to the Regge- Wheeler potential ^ 
of a Schwarzschild black hole. Thus, the axial potential barrier generated by the curvature of 
the star depends on how the energy-density and the pressure are distributed inside the star in 
the equilibrium configuration, and therefore it depends on the equation of state of matter inside 
the star As an example, in Figure |2] we show the € - 2 potential barrier for an ideal, constant 
density star with /?/M = 2.8 (left panel) and RjM - 2.4 (right panel). If we compare the potential 
shown in Figure |2] with the Regge-Wheeler potential of a Schwarzschild black hole shown in 
Figure [T] we notice an important difference. The Schwarzschild potential vanishes at the black 
hole horizon, and has a maximum at r,„ax ~ 3M, whereas the potential barrier of a perturbed 
star tends to infinity at r = 0. Thus, for a Schwarzschild black hole waves are scattered by a 
one-dimensional potential barrier, whereas in the case of a star they are scattered by a central 
potential. 



Since the axial perturbations do not excite any motion in the fluid, for a long time they have 
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Figure 2. The £ = 2 potential barrier for a constant density star of mass M = 1.4 Mq and R/M = 2.8 (left 
panel) and R/M = 2.4 (right panel). The horizontal, dashed line in the right panel corresponds to the value 
of frequency a/, which corresponds to a solution regular at r = and behaving as a pure outgoing wave at 
radial infinity, i.e. to a quasi-normal mode. 

been considered as trivial. But this is not true if we adopt the scattering approach: the absence 
of fluid motion simply means that the incident axial wave experiences a potential scattering, and 
this scattering can, in some extreme conditions, be resonant. Indeed, if we look for solutions that 
are regular at r = and behave as pure outgoing waves at infinity, we find modes which do not 
exist in Newtonian theory; if the star is extremely compact, the potential in the interior is a well, 
and if this well is deep enough there can exist one or more more slowly damped quasi-normal 
modes, or .s-modes (Chandrasekhar & Ferrari 1991c). For example, if the mass of the star is, say, 
M = 1.4 Mq and R/M - 2.4, i.e. the stellar compactness is M//? - 0.42, as shown in Figure |2] 
the well inside the star is deep enough to allow one quasi-normal mode. The number of s-modes 
increases with the depth of the well, which corresponds to a larger stellar compactness. However, 
it should be mentioned that neutron stars are not expected to have such a large compactness, 
unless one invokes some exotic equation of state. The .s-modes are also named trapped modes 
because, due to the slow damping, they are effectively trapped by the potential barrier, and not 
much radiation can leak out of the star when these modes are excited. Axial modes on a second 
branch are named w-modes and are highly damped (Kokkotas 1994). The w-mode frequency 
also depends on the stellar compactness, as we shall show in Section [TTI Therefore they carry 
interesting information on the internal structure of the star. 

It should be stressed that the axial modes do not have a Newtonian counterpart. 

Our approach to the polar perturbations, which couple the perturbations of the gravitational 
field to those of the metric, is different from the Newtonian approach briefly described in the in- 
troduction. Rather than focusing on the fluid behaviour, we focus on the variables which describe 
the spacetime perturbations, assuming that, as in the case of black holes, they are excited by 
the incidence of polar gravitational waves belonging to a particular angular harmonic. A careful 
scrutiny of the structure of the polar equations shows that it is possible to decouple the equations 
describing the metric from those describing the fluid perturbations. This decoupling allows us to 
solve the equations for the spacetime perturbations with no reference to the motion that can be in- 
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duced in the fluid, and fliis is possible in general. Once flie solution for the metric perturbations is 
found, the fluid variables can be determined in terms of them by simple algebraic relations with- 
out further ado (Chandrasekhar & Ferrari 1990b). The final set of equations to solve is described 
in Section 2.1. 



2.1 The equations for the polar perturbations 



Assuming that the metric which describes the unperturbed star has the form 



■ e'^dtf-e^dip- 



{drf 



«2» 



{def. 



(13) 



the functions that describe the polar perturbations, expanded in spherical tensor harmonics and 
Fourier-expanded are 



6v = Nt{r)P ({cofi 9)e'"' 

6p3 = lTeir)Pt + V((r)P(,e,ey' 
6p = n^OPKcos e)e'"' 
6e = Ee{r)Pe{cos 9)e'"' 



6xlJ = [Te{r)P{ + Ve(r)P e,e cot 6]e'"' , 
lie + p)e'^^'^^,.{r, 0)e''" = Uiir)P(e'"' 
2(e + /,)e''+'"ffl(r,0)e'"' = WMPt^" 



(14) 



where (cos 6) are Legendre's polynomials, to is the frequency, 6p and 5e are perturbations of 
the pressure and of the energy density, and ^r, are the relevant components of the Lagrangian 
displacement of the generic fluid element. Note that (A^, L, T, V) and (11, E, U, W) are, respec- 
tively, the radial part of the metric and of the fluid perturbations. After separating the variables 
the relevant Einstein's equations for the metric functions become 

(r^Gi),,. = nvANr - U) + ^(^2/^2 _ + + r{v^, - p2,r)Xe,r + o?e'^^^-''^rXt, 
-v,rNe,r = -Ge + v,AX(,r + y/Nc - Le)] + ^.ie^"^ - \)(Ne - rX,,, - r^Gt) 
-e^f^ie + p)Nt + i^^2g2te-v) + + + L[^xe,r + (2n + l)Xi]} , 

Lf,,(l -D)+L( [(^ - y - (i + V ,) d] + Xt,r +Xt{\- V,r) + DNt,r + 

Nt {Dvr -T-P) + {-r+ Ev,r) [n^ - L( + '^Gi + \ (rX,,, + X()\ = 0, 



(15) 



where 



Q 



P I enlrops-consl ' 



(£+£) 
JV ' 



B 



2(e+p) 



(fr - QV,y)^ 



life (g+p) 



{A-- 
1 -- 

E^D{Q-\)-Q, 

2(A+B) 2oj^e-^''(e+p)+e-^2v^(e_r-Qp^) ' 

and V( and T( have been replaced by Xe and Gf defined as 

tXf = nVt 

Gc = vA^^Xt - Tfl,, + ^(e^/'^ - l)[n(N( + Tt) + Ne] 
[ V-f{N( + Le) - e2/'2(e + p)Ne + \a?e^'^--^\L, -Te+ ^Xd 



(16) 



(17) 
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These equations are valid in general, also for non-barotropic equations of state. It should be 
stressed that equations ( fTSl l govern the variables (X,G,N,L) which are metric perturbations; 
however, since the motion of the fluid is excited by the polar perturbation, we may want to 
determine the fluid variables, (O, E, U, W); they can be obtained in terms of the metric functions 
using the following algebraic relations 

Tie = -^co^e-^^'Wi - (e + p)Nr, Eg = QUe + ^ (c, - Qp.rWc, 
2 2(e + p) 

[(a?e-^'W(),r + (Q+ IXricJ^e-^^Wt) + 2(e, - Qp,r)Nr](€ + p) 

Up — 

Ke-2''(e + p) + e-2/^^v,(e, - gp,,)] 

Outside the star the fluid variables vanish, and the polar equations reduce to the wave equation 
© with the Zerilli potential ©. 

As discussed above, for the axial perturbations, the frequencies of the quasi-normal modes 
were found by solving a problem of scattering by a central potential; for the polar perturbations it 
is not so simple, because a Schroedinger equation holds only in the exterior of the star, whereas 
a higher order system must be solved in the interior. It is still a scattering problem, but of a more 
complex nature since the incident polar gravitational waves, which excite the perturbations, drive 
the fluid pulsations, which in turn emit the scattered component of the wave. This approach was 
very fruitful in many respects. First of all, given the equilibrium configuration for any assigned 
equation of state, it was very easy to evaluate the QNM-frequency by integrating the equations 
for the metric perturbations inside and outside the star, looking for the solutions which, being 
regular at r = 0, behave as pure outgoing waves at infinity. Furthermore, we generalized the 
perturbed equations to slowly rotating stars, and derived the equations which describe how the 
axial perturbations couple to the polar (Chandrasekhar & Ferrari 1991b). 



2.2 Perturbed equations for a slowly rotating star 

Very briefly, the coupling mechanism is the following. Let be the axial radial function, 
solution of equation ( fTOl i. which describes the perturbation of a non-rotating star; let e(0)Z'" be 
the perturbation to first order in the star angular velocity O. The axial perturbation is the sum of 
the two: 

As Zj?", the function satisfies the wave equation (fTOt with the same potential (fTSl i. but with a 
forcing term: 

(=2 [ * j e=2 

3 

where p = cos-d- and C^^^il^) ^i"^ Gegenbauer polynomials. The source term 5^ is 
5° = Z!T_,[(2W° -I- A^J? + 51}} + 2nV^P(,^ + 2pV°Pe,^,^] + 2mW°(Q - IXrPe,/, 
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it is a combination of the functions which describe the polar perturbations on the non-rotating 
star, found by solving the equations given in Section IZTI It should be stressed that the coupling 
function m is the function responsible for the Lense-Thirring effect. Thus a rotating star exerts 
a dragging not only of the bodies, but also of the waves, and consequently an incoming polar 
gravitational wave can convert, through the fluid oscillations it excites, some of its energy into 
outgoing axial waves. This is a purely relativistic effect, and it is due to the dragging of inertial 
frames. It is interesting to note that the coupling between axial and polar perturbations satisfies 
rules that are similar to those known in the theory of atomic transitions: a Laporte rule and a 
selection rule, according to which the polar modes belonging to even ( can couple only with the 
axial modes belonging to odd t, and conversely, and that it must be 

I - m+\, or I - in - 1. 

Furthermore, the coupling satisfies a propensity rule (Fano 1985); the transition ^ — > ^ + 1 is 
strongly favoured over the transition £ { - 1. 

At the time Chandrasekhar and I wrote the series of papers on stellar perturbations, there 
was a growing interest in the subject, also motivated by the fact that the construction of ground 
based interferometric detectors, LIGO in the US and Virgo in Italy, had just started. Many studies 
addressed the problem of finding the frequencies of the QNMs, to establish what kind of infor- 
mation they carry on the internal structure of the emitting source. The collective effort developed 
using essentially two different perturbative approaches: one in the frequency domain, as for the 
theory developed by Thorne and collaborators or by Chandrasekhar and myself, another in the 
time domain. The time domain approach basically consists in separating the equations of stellar 
perturbations as usual in terms of spherical harmonics, and in solving the resulting equations in 
terms of two independent variables, radial distance and time. The equations are excited using 
some numerical input, like for instance a Gaussian impulse, and then the QNM frequencies are 
found by looking at the peaks of the Fourier transform of the signal obtained by evolving the 
time-dependent equations numerically. A disadvantage of this evolution scheme is that one can- 
not get the complete spectrum of the QNMs either for a star or for a black hole. The reason is 
that, although any perturbation is the sum of the harmonics involved, in practice only a few of 
them can be clearly identified; thus, to find some more modes one has to proceed empirically by 
changing the initial conditions. However, the evolution of the time dependent equations is, still 
today, the only viable perturbative method to find the QNM frequencies and waveforms emitted 
by rapidly rotating relativistic stars. To describe the problems which emerge when dealing with 
the perturbations of a rapidly rotating star is beyond the scope of this paper; I will discuss some 
related issues in the concluding remarks. 



3. Neutron star oscillations 



We shall now show how the theory of perturbations of non-rotating stars can be applied to gain 
some insight into the internal structure of the emitting source. Different classes of modes probe 
different aspects of the physics of neutron stars. For instance the fundamental mode (/-mode). 
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which has been shown to be the most efficient GW emitter by most numerical simulations, de- 
pends on the average density, the pressure modes (p-modes) probe the sound speed throughout the 
star, the gravity modes (g-modes) are associated with thermal/composition gradients and the w- 
modes are spacetime oscillations. Furthermore, crustal modes, superfluid modes, magnetic field 
modes can, if present, add to the complexity of stellar dynamics. The sensitivity of ground based 
gravitational detectors has steadily improved over the years in a broad frequency window; the 
advanced version of LIGO and Virgo, and especially third generation detectors like ET, promise 
to be powerful instruments to detect signals emitted by oscillating stars. The frequencies of quasi 
normal modes are encoded in these signals; therefore, as the Sun oscillation frequencies are used 
in helioseismology to probe its internal structure, we hope that in the future it will be possible to 
use gravitational waves to probe the physics of neutron stars. One of the issues which is interest- 
ing to address concerns the equation of state of matter in a neutron star core, which is actually 
unknown. This problem is of particular interest, because the energies prevaihng in the inner core 
of a neutron star are much larger than those accessible to high energy experiments on Earth. In 
the core, densities typically exceed the equilibrium density of nuclear matter, po - 2.67 x 10'^ 
g/cm^; at these densities neutrons cannot be treated as non-interacting particles, and the main 
contribution to pressure, which comes from neutrons, cannot be derived only from Pauli's exclu- 
sion principle. Indeed, with only this contribution, we would find that the maximum mass of a 
neutron star is 0.7 Mq which, as observations show, is far too low. This clearly shows that NS 
equilibrium requires a pressure other than the degeneracy pressure, the origin of which has to 
be traced back to the nature of hadronic interactions. Due to the complexity of the fundamental 
theory of strong interactions, the equations of state appropriate to describe a NS core have been 
obtained within models, which are constrained, as much as possible, by empirical data. They 
are derived within two main, diff'erent approaches: the nonrelativistic nuclear many-body theory, 
NMBT, and the relativistic mean field theory, RMFT. In NMBT, nuclear matter is viewed as a 
collection of pointUke protons and neutrons, whose dynamics is described by the nonrelativistic 
Hamiltonian: 

^=Z^+Z^o + Z^o*. (19) 

i j>i k> j>i 

where m and pi denote the nucleon mass and momentum, respectively, whereas v,y and Vijk de- 
scribe two- and three-nucleon interactions. These potentials are obtained from fits of existing 
scattering data (Wiringa, Stoks & Schiavilla 1995), (Pudliner et al. 1995). The ground state en- 
ergy is calculated using either variational techniques or G-matrix perturbation theory. The RMFT 
is based on the formalism of relativistic quantum field theory, nucleons are described as Dirac 
particles interacting through meson exchange. In the simplest implementation of this approach 
the dynamics is modeled in terms of a scalar and a vector field (Walecka 1974). The equations 
of motion are solved in the mean field approximation, i.e. replacing the meson fields with their 
vacuum expectation values, and the parameters of the Lagrangian density, i.e. the meson masses 
and coupling constants, can be determined by fitting the empirical properties of nuclear matter, 
i.e. binding energy, equilibrium density and compressibility. Both NMBT and RMFT can be 
generalized to take into account the appearance of hyperons. In the following we shall consider 
some EOS representative of the two approaches, which have been used in the literature. 
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It should be stressed that different ways of modehng hadronic interactions affect the pulsation 
properties of a star, which we are going to discuss. 



3.1 The axial and polar w-modes 

As shown in Section |2] the axial perturbations are described by a Schroedinger-like equation 
with a central potential barrier which depends on the energy and pressure distribution in the un- 
perturbed star, i.e. on the equation of state. The slowly damped modes are not expected to be 
associated with significant gravitational wave emission, because they are effectively trapped by 
the potential barrier; in addition they appear if the star has a compactness close to the static 
Schwarzschild limit, which establishes that constant density star solutions of Einstein's equations 
exists only for M/R < 4/9 ^ 0.44 . Conversely, the w-modes, which are highly damped, ex- 
ist also for stars with ordinary compactness. They have been shown to exist also for the polar 
perturbations and in that case they are coupled to negligible fluid motion. In Figure |3] we com- 
pare the frequencies of the lower axial vv-modes computed in (Benhar, Berti & Ferrari 1999) 
with those of the lower polar vv-modes computed in (Andersson & Kokkotas 1988) for several 
EOSs. The main features of different EOS are, very briefly, the following. EOS A (Pandhari- 
pande 1971a) is pure neutron matter, with dynamics governed by a nonrelativistic Hamiltonian 
containing a semi-phenomenological interaction potential. It is obtained using NMBT EOS B 
(Pandharipande 1971b) is a generalization of EOS A, including protons, electrons and muons in 
j0-equilibrium, as well as heavier baryons (hyperons and nucleon resonances) at sufficiently high 
densities (NMBT). EOS WFF (Wiringa, Fiks & Fabrocini 1988) is a mixture of neutrons, protons, 
electrons and muons in y6-equilibrium. The Hamiltonian includes two- and three-body interaction 
potentials. The ground state energy is computed using a more sophisticated and accurate many- 
body technique (NMBT). In EOS L (Pandharipande & Smith 1975) neutrons interact through 
exchange of mesons (Lij,p, cr). The exchange of heavy particles (w,p) is described in terms of 
nonrelativistic potentials, the effect of cr-meson is described using relativistic field theory and the 
mean-field approximation. 

From Figure |3] we see that for each selected EOS the frequency of the polar w-modes is a 
rather steeply decreasing function of the stellar compactness M/R, whereas for the axial modes 
the dependence of v,,^, on the compactness is weak, and ranges within intervals that are separated 
for each EOS. This means that if an axial gravitational wave emitted by a star at a given frequency 
could be detected, we would be able to identify the equation of state prevailing in the star's interior 
even without knowing its mass and radius. Hence, the detection of axial gravitational waves 
would allow us to constrain the EOS models, with regard to both the composition of neutron star 
matter and the description of the hadronic interactions. Until very recently, the common belief 
was that w- modes are unlikely to be excited in astrophysical processes. However, it has been 
shown that they are excited in the collapse of a neutron star to a black hole, just before the black 
hole forms (Baiotti et al. 2005). Unfortunately the typical frequencies of these modes (of the 
order of several kHz) are higher than the frequency region where the actual gravitational wave 
detectors are sensitive. 
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Figure 3. The frequency of the first polar (dashed line) and axial (continuous line) w-modes are plotted as 
a function of the star compactness for the EOSs A, B, WFF, L. 



3.2 Polar Quasi Normal Modes 



The polar metric perturbations are physically coupled to the fluid perturbations. As shown in Sec- 
tion 12.11 the frequencies of the polar QNMs can be computed by solving a system of equations 
involving only the metric perturbations; however, they carry a strong imprint of the internal com- 
position of the star, which is present in equations (flST l through the pressure and energy density 
profiles in the unperturbed star, which appear as coefficients of the differential equations. Ac- 
cording to a scheme introduced by Cowling in Newtonian gravity (Cowling 1942), polar modes 
can be classified on the basis of the restoring force which prevails when the generic fluid element 
is displaced from the equilibrium position: for ^-modes, or gravity modes, the restoring force is 
due to buoyancy, for p-modes it is due to pressure gradients. The mode frequencies are ordered 
as follows 

■ Mg^^ < .. < (Ajg^ < (Jjf < LOp^ < .. < OJp,^-. 

and are separated by the frequency of the fundamental mode (/-mode), which has an intermediate 
character between g- and p- modes. As discussed in Section lTTl general relativity predicts also 
the existence of polar w-modes, that are very weakly coupled to fluid motion and are similar to 
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the axial w-modes (Kokkotas & Schutz 1992). The frequencies of axial and polar w-modes are 
typically higher than those of the fluid modes g, f and p. 

If we are mainly interested in gravitational wave emission, the most interesting mode is the 
/-mode. For mature neutron stars, its frequency is in the range 1-3 kHz, which is in the 
bandwidth of ground based detectors Virgo and LIGO (although not in the region where they 
are most sensitive); the damping times are of the order of a few tenths of seconds, therefore the 
excitation of the /-mode would appear in the Fourier transform of a gravitational wave signal as 
a sharp peak and could, in principle, be extracted from the detector noise by an appropriate data 
analysis. Moreover, the fundamental mode could be excited in several astrophysical processes, 
for instance in the aftermath of a gravitational collapse, in a glitch, or due to matter accretion onto 
the star. For this reason, since the early years of the theory of stellar perturbations, the interest 
of scientists working in this field has initially been focussed on the determination of the /-mode 
frequencies. After the work of Lindblom & Detweiler in 1983 and of Cutler & Lindblom in 1987, 
who respectively computed the t - 2 and H > 2 /-mode eigenfrequencies for the EOSs available 
at that time, more recently this work has been updated, and extended to other modes, by Anderson 
& Kokkotas (1998) and Benhar, Ferrari & Gualtieri (2004). In particular, in these two papers the 
/-mode frequency Vf and the corresponding damping time tj- have been computed to establish 
whether V/ scales with the average density of the star, as it does in Newtonian gravity, and whether 
there also exists a scaling law for Ty . The sets of EOSs used in the two works are not identical, 
because the papers were written six years apart, although some EOSs appear in both (see the two 
papers for details). The work done by Benhar, Ferrari & Gualtieri (2004) also includes examples 
of hybrid stars, namely neutron stars with a core composed of quarks. In Anderson & Kokkotas 
(1998) Vf and t/ have been fitted by a linear function of the average density of the star (M/R^y^^, 
and of its compactness M/R, as follows. 



0.78 -H 1.635 



1 



22.85 - 16.65 



(20) 



where M - M/1.4 Mq and R = R/{10 km). Here and in the following formulae v/ is expressed 
in kHz and t/ in s. The fits for vy and t/ obtained by Benhar, Ferrari & Gualtieri (2004) using 
the new set of EOSs are 



a = 0.79 ± 0.09 (in kHz), = 33 ± 2 (in km). 



(21) 



and 



1 



cM^ r (M) 
— T- a + b\ 
R^ 



a = [8.7 ±0.2] ■ 10" 



b = -0.271 ±0.009. 



(22) 



In equations (|2TI) and (l22l i mass and radius are in km (i.e. mass is multiplied by G/c^) and 
c = 3 ■ 10^' km/s. The data for the different EOSs used by Benhar, Ferrari & Gualtieri (2004), 
and the fits given in equations (I20]i- (l22l i are shown in Figure |4] vy is plotted in the upper panel 
as a function of the average density; the fit (l20l i is shown as a black dashed line labelled 'AK fit'. 
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whereas the new fit is indicated as a red continuous Hne labelled 'NS fit' . The NS fit is lower 
by about 100 Hz than the AK fit, showing that the new EOSs are, on average, less compressible 
than the old ones. The quantity {R^ j cM^) I t j given in equation ( l22l i is plotted in the lower panel 
of Figure H] versus the stellar compactness MjR. In this case the AK fit for r/ ( l20t . and the NS 
fit (I22I 1 are nearly coincident. For comparison, in both panels of Figure |4] we show the frequency 
and the damping time of the /-mode of a population of strange stars, namely stars entirely made 
of up, down and strange quarks, modeled using the MIT Bag-model, spanning the allowed range 
of parameters, which are the Bag constant, the coupling constant as and the quark masses (see 
Benhar et al. 2007 for details). The parameters of the fits for strange stars are 

for Vf a = -[0.8 + 0.08] ■ 10"^ , b^A6 + 0.2, (23) 

and 

for r/ fl = [4.7 ± 5 ■ 10"^] ■ 10"^ ^7 = -0.12 ± 3 ■ 10"^ (24) 

In Figure |4] the fits for strange stars are labelled as 'SS fit'. It is interesting to note that the SS 
fits are quite different from those appropriate for neutron stars (AK- and NS-fits). First of all the 
errors on the parameters are much smaller, which indicate that the linear behaviour is followed 
by these stars, both for vj and for Tf, irrespective of the values of the parameters of the model. 
Moreover, the difference between the fits is much larger for lower values of the average density. 

The empirical relations given in eqs. (I20t- (l24t could be used to constrain the values of the 
star mass and radius, were the values of vy and tj identified in a detected gravitational signal. 
The stellar parameters would be further constrained if other modes are excited and detected and, 
knowing them, we would gain information on the equations of state of matter in the neutron star 
core, whose uncertainty is due, as explained earlier, to our ignorance of hadronic interactions. 
Furthermore, if the neutron star mass is known, as it may be if the star is in a binary system, 
the detection of a signal emitted by the star oscillating in the /-mode may provide some further 
interesting information (Benhar et al. 2007). In Figure |5] we plot v/ as a function of the stellar 
mass, for neutron/hybrid stars and for strange stars modeled using the MIT bag model. Note that 
1.8 Mq is the maximum mass above which no stable strange star can exist. We see that there is 
a small range of frequency where neutron/hybrid stars are indistinguishable from strange stars; 
conversely, there is a large frequency region where only strange stars can emit. Moreover, strange 
stars cannot emit gravitational waves with vj < 1.7 kHz, for any value of the mass in the range 
we consider For instance, if the stellar mass is M = 1.4 M©, a signal with vj > 2 kHz would 
belong to a strange star. Figure |5] also shows that, even if we do not know the mass of the star 
(as it is often the case for isolated pulsars), if v/ > 2.2 kHz, apart from a very narrow region of 
masses where stars with hyperons would emit (EOS BBS 1 and G240), we can reasonably rule out 
that the signal is emitted by a neutron star In addition, it is possible to show that, since v/ is an 
increasing function of the Bag constant B, if a signal emitted by an oscillating strange star were 
detected, it would be possible to set constraints on B much more stringent than those provided by 
the available experimental data (Benhar et al. 2007). 

In conclusion, the QNM frequencies can be used to gain direct information on the equation 
of state of matter in a neutron star core. 
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Figure 4. The frequency of the fundamental mode is plotted in the upper panel as a function of the square 
root of the average density for the different EOSs considered by Benhar, Ferrari & Gualtieri (2004). We 
also plot the fit given by Anderson & Kokkotas (1998) plotted as AK-fit and our fit (NS-fit). The NS-fit is 
systematically lower (about 100 Hz) than the AK-fit. The damping time of the fundamental mode is plotted 
in the lower panel as a function of the compactness MjR. The AK-fit and our fit, plotted respectively as a 
dashed and continuous line, do not show significant differences. 



18 



V. Ferrari 







APR2 




Neutron stars 


APRB200 -■- 






APRB120 — ^ ' 




+ 


BBS1 — •— 




y \ Strange stars 


G240 - 









0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 

M/Mq 

Figure 5. The frequency of the fundamental mode is plotted as a function of the mass of the star, for 
neutron/hybrid stars (continuous lines) and for strange stars modeled using the MIT bag model, spanning 
the set of parameters indicated in the range allowed by high energy experiments (dashed region). 



The crucial question now is: do we have a chance to detect a signal emitted by a star oscillat- 
ing in a polar quasi-normal mode? Detection chances depend on how much energy is channeled 
into the pulsating mode, which is unknown, and on whether the mode frequency is in the detector 
bandwidth. The signal emitted by a star pulsating in a given mode of frequency v and damping 
time T, has the form of a damped sinusoid 

h(t) = Jle-^'-'"^'^ sin[27rv(f - fo)] for t > to, (25) 

where fo is the arrival time of the signal at the detector (and h{t) - for f < fo). The wave 
amplitude ^ can be expressed in terms of the energy radiated in the oscillations. 



where AEq = AEcw/MqC^. This quantity is unknown. Therefore to assess the detectability of a 
signal we can only evaluate how much energy should be emitted in a given mode, in order for the 
signal to be detected by a given detector with an assigned signal to noise ratio {S /N) 

(If.^^ (27) 
\n) 1+422 25/ ^ ' 

In this equation Q = ttvt is the quality factor and S„ is the detector spectral noise density. As- 
suming as a bench-mark for AEq the energy involved in a typical pulsar glitch, in which case a 
mature neutron star might radiate an energy of the order of AE^w = 10^^^ MqC^ , and assuming 
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V ~ 1500 Hz, T ~ OA s, d = 1 kpc, we find J?[ » 5 x 10~ . Such a signal is too weak to be 
seen by actual detectors, therefore we conclude that 3rd generation detectors are needed to detect 
signals from old neutron stars. More promising are the oscillations of newly-born neutron stars; 
indeed, since a NS forms as a consequence of a violent, and generally non-symmetric event - the 
gravitational collapse - a fraction of its large mechanical energy may go into non-radial oscilla- 
tions and would be radiated in gravitational waves. Thus, during the first few seconds of the NS 
hfe more energy could be stored in the pulsation modes than when the star is cold and old. In 
addition, during this time the star is less dense than at the end of the evolution; consequently, the 
frequencies of the modes which depend on the stellar compactness (as for instance the /-mode) 
are lower and therefore span a frequency range were the detectors are more sensitive (Ferrari, 
Miniutti & Pons 2003). For instance, if we assume that an energy AEgw = 1.6 • 10"^ MqC^ is 
stored in the /-mode of a neutron star just formed in the Galaxy, the emitted signal would be 
detectable with a signal to noise ratio S/N by advanced Virgo/LIGO, and with S/N = 2.7 by 
Virgo-i-/LIGO, the upgraded configurations now in operation. 

4. Stellar perturbations and magnetar oscillations 

Magnetars are neutron stars whose magnetic field is, according to current models, as large as 
lO'^ G (Thompson & Duncan 1993, 2001). During the last three decades some very interesting 
astrophysical events have been observed which are connected to magnetar activity and stellar 
pulsations. They involve Soft Gamma Repeaters (SGRs), which are thought to be magnetars; 
these sources occasionally release bursts of huge amount of energy (L ^ lO''^ - 10^^ ergs/s), and 
these giant flares are thought of being generated from large-scale rearrangements of the inner 
field, or catastrophic instabilities in the magnetosphere (Thompson & Duncan 2001; Lyutikov 
2003). Up to now, fliree of these events have been detected: SGR 05026-66 in 1979, SGR 
1900-Hl4in 1998 and SGR 1806-20 in 2004. In twooffliem (SGR 1900+14 and SGR 1806-20), 
a tail lasting several hundred seconds has been observed, and a detailed study of this part of the 
spectrum has revealed the presence of quasi-periodic oscillations (QPOs) with frequencies 

18, 26, 30, 92, 150, 625 and 1840 Hz 

for SGR 1806-20 (Watts & Strohmayer 2006), and 

28, 53, 84 and 155 Hz 

for SGR 1900-1-14 (Strohmayer & Watts 2006). The discovery of these oscillations stimulated an 

interesting and lively debate (still ongoing) among groups working on stellar perturbations, about 
the physical origin of these sequences. Of course in order to study these oscillations the magnetic 
field and its dynamics have to be included in the picture, while rotation plays a less important role, 
since observed magnetars are all very slowly rotating. The problem presents extreme complexity 
both at conceptual and at computational levels; therefore it is usually approached using simplify- 
ing assumptions and/or approximations. Some of the studies try to explain the observed modes 
in terms of torsional oscillations of the crust (Samuelsson & Andersson 2007; Sotani, Kokkotas 
& Stergioulas 2007), others attribute the observed spectra to global magneto-elastic oscillations 
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(Glampedakis, Samuelsson & Andersson 2006), still others investigate the interaction between 
the torsional oscillations of the magnetar crust and a continuum of magnetohydrodynamic modes 
(the Alfven continuum) in the fluid core (Levin 2007; Sotani, Kokkotas & Stergioulas 2008; Co- 
laiuda, Beyer & Kokkotas 2009; Cerda-Duran, Stergioulas & Font 2009; Colaiuda & Kokkotas 
2010; Gabler et al. 2011) using different approaches and approximations. In particular, in (Co- 
laiuda & Kokkotas 2010) the torsional oscillations of a magnetar have been studied in a general 
relativistic framework, perturbing Einstein's equations in the Cowling approximation, i.e. ne- 
glecting gravitational field perturbations. By this approach the crust-core coupling due to the 
strong magnetic field has been shown to be able to explain the origin of the observed frequen- 
cies, at least for SGR 1806-20, if a suitable stellar model is considered. With this identification, 
constraints on the mass and radius of the star, and consequently on the EOS in the core, can be 
set; estimates of the crust thickness and of the value of the magnetic field at the pole can also be 
inferred. 

Thus, the theory of stellar perturbations has been generalized to magnetized stars, although 
for now only with considerable restriction, since only torsional oscillations have been considered 
(i.e. axial perturbations) and only in the Cowling approximation. Nevertheless, it already pro- 
vides very interesting information on the dynamics of these stars and allows us to confront the 
predictions with astronomical observations. 



5. Concluding remarks 



I would like to conclude this review by mentioning the fact that the theory of perturbations of 
rotating stars has not been developed to the same extent as the theory of non-rotating stars. The 
main reason is that the mathematical tools appropriate for a successful variable separation has not 
been found yet. When the perturbations of a non-rotating black hole are studied, separation of 
variables is achieved by expanding all tensors in tensorial spherical harmonics. In the case of Kerr 
perturbations, namely of perturbations of an axisymmetric, Petrov type D background, the same 
result is obtained by expanding the Newman-Penrose quantities in oblate spheroidal harmonics. 
When perturbing a rotating star, i.e. an axisymmetric solution of Einstein-i-hydro equations, an 
expansion in terms of tensorial spherical harmonics leads, as we have seen in the case of slow 
rotation in Section l272l to a coupling between polar and axial perturbations. If rotation is not slow, 
the number of couplings to be considered increases to such an extent that the problem becomes 
untreatable, both from a theoretical and from a computational point of view. One may argue that, 
since the background of a rapidly rotating star is not spherically symmetric, tensorial spherical 
harmonics are unappropriate, and this is certainly true. However, even if we try to use the oblate 
spheroidal harmonics and the Newmann Penrose formalism we fail: the coupling between the 
metric and the fluid makes the separation impossible, at least in terms of these harmonics (unlike 
the Kerr metric, the metric describing a star is not of Petrov type D). For this reason, perturba- 
tions of rotating stars have been studied either in the slow rotation regime, or using the Cowling 
approximation, which neglects spacetime perturbations, or using other simplifying assumptions. 
For instance, as far as the mode calculation is concerned, the Cowling approximation allows de- 
termination with reasonable accuracy the frequency of the higher order p-modes, of the g-modes 
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and of the inertial modes, like the r-modes, thus allowing us to gain information on the onset of 
related instabilities. Conversely, the determination of the /-mode frequency, which is so impor- 
tant from the point of view of gravitational wave emission, is not very precise, leading to errors 
as large as ~2Q%. 

However, it should be mentioned that non-linear simulations of rotating stars are producing 
very interesting results; for instance in a recent paper Stergioulas and collaborators (Zink et al. 
2010) have been able to follow the frequency of the non-axisymmetric fundamental mode of a 
sequence of rotating stars with increasing angular velocity, up to the onset of the CFS instabil- 
ity, making also very optimistic estimates of the amount of gravitational radiation which could 
be emitted in the process. To describe matter in the neutron star they use a simple model (a 
polytropic equation of state and uniform rotation); however, their result indicates that numerical 
relativity is making giant steps in this field. Thus, supercomputers are making accessible very 
complex problems, which only ten years ago one would not have dreamed of solving; however, 
perturbation theory still remains a very powerful tool to investigate many physical problems and 
it should be used in parallel with the numerical work to gain a deeper insight into the physics of 
stellar oscillations. 
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